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Abstract. The internal dynamics of strongly interacting systems and that of bio-polymers 
such as proteins display several important analogies, despite the huge difference in their 
characteristic energy and length scales. For example, in all such systems, collective excitations, 
cooperative transitions and phase transitions emerge as the result of the interplay of strong 
correlations with quantum or thermal fluctuations. In view of such an observation, some 
theoretical methods initially developed in the context of theoretical nuclear physics have been 
adapted to investigate the dynamics of biomolecules. In this talk, we review some of our 
recent studies performed along this direction. In particular, we discuss how the path integral 
formulation of the molecular dynamics allows to overcome some of the long-standing problems 
and limitations which emerge when simulating the protein folding dynamics at the atomistic 
level of detail. 



1. Introduction 

Proteins are macromolecules involved in virtually all the biochemical processes which take 
place inside the cell, such as e.g. catalysis, charge transport, signal transduction. Clearly, 
understanding the physical principles which drive the internal dynamics of individual proteins 
and shape the protein-protein interaction would have countless implications in molecular biology, 
drug design and nano-biotechnology. In particular, a central open challenge at the interface 
of physics, biochemistry and molecular biology concerns the prediction of the protein folding 
pathways i.e. the sequence of conformational transitions which these molecules perform in order 
to reach their stable and biologically active configuration (see e.g. Fig. [TJ . 

From a theoretical perspective, the problem of determining the time evolution of proteins has 
been mostly investigated by means of Molecular Dynamics (MD) simulations, i.e. through the 
direct integration of the classical equations of motion of all the atoms in the system. From the 
conceptual point of view, such an approach should provide a reasonable approximation, since 
quantum corrections are expected to be quite small. On the other hand, from the practical 
point of view, the applicability of classical MD simulations is strongly limited by their intrinsic 
computational complexity. Indeed, even using computational resources comparable with those 
employed in the largest contemporary lattice QCD simulations, it is possible to follow the 
dynamics over time intervals which do not exceed a few hundreds ns [lj. Unfortunately, the 
dynamics which drives the protein phenomenology occurs at time scales which are several orders 
of magnitude larger, typically ranging from few ms to several seconds. 
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Figure 1. Schematic representation of a protein folding pathway. This trajectory corresponds 
to a biological instanton and was obtained using the Dominant Reaction Pathways approach 
(see discussion in section [4j 

The decoupling of the microscopic time scales from the time scales at which conformational 
reactions occur is a common feature in many bio-molecular systems. It is due to the fact 
that most conformational transitions are cooperative and involve the overcoming of free-energy 
barriers. As a result, in a MD simulation, most computational time is wasted to describe thermal 
oscillations inside the initial meta-stable state. 

Given the huge gap between the time scales at which the relevant physics takes place and 
the time intervals which can be covered by MD simulations, it is natural to expect that a real 
breakthrough in our comprehension of biological matter at the fundamental level is more likely 
to come from new ideas and new theoretical approaches, rather than from the next several 
generations of supercomputers. 

In the next sections, we discuss how functional integral approaches originally developed 
in the context of nuclear theory can provide an extremely powerful framework to investigate 
the dynamics of proteins as well as other complex molecular systems. Before introducing the 
mathematical details of this theory, it is instructive to discuss some basic aspects of protein 
structure and phenomenology and to highlight several analogies with the dynamics of strongly 
interacting systems. 

2. QCD and Protein Dynamics: So Far Yet So Close 

From the biochemical point of view, proteins are poly-peptide chains, i.e. weakly branched 
polymers formed by 20 different types of monomers called amino- acids. Typical globular proteins 
consist of about 100 amino- acids, but some chains can contain as many as ~ 30, 000 of such 
monomers. 

From the physical point of view [2J, proteins are very different from essentially all other 
known polymers: indeed, while standard hetero-polymers are prototypes of disordered mesoscale 
systems, proteins are characterized by their unique "ability" to spontaneously and reversibly fold 
into a well-defined stable configuration (denominated the native state), in which they perform 
their biological functions (see e.g. Fig. [2]). Such a remarkable feature is the result of the 
evolutionary pressure, which has selected only a very specific subset in the space of all possible 



Figure 2. Ball-and-stick representation of the native state of a globular protein. The tube 
represents the back-bone of the chain. 

amino- acid sequences. 

The native state of a protein is characterized by an extremely low entropy, since it is defined 
by the small thermal oscillations around a single conformation. In this sense, proteins may be 
regarded as non-symmetric mesoscale crystals. If the temperature of the water surrounding the 
molecule is increased, or if a denaturant agent (typically urea) is added to the solvent, then 
proteins can reach a much more entropic phase called the molten globule. In this state, the 
chain is still collapsed, but is free to explore the different compact configurations, in analogy 
with a liquid drop. Finally, at even higher temperatures or denaturant concentrations, proteins 
swell and start to behave as flexible random coils (see right panel of Fig. [3]). 

Interestingly, the phase diagram of proteins in solution displays several analogies with that 
of strongly interacting matter (see left panel of Fig. [3]). Indeed, the low-temperature and 
low-density phase of QCD is characterized by a very low entropy, since only colorless hadrons 
are present in the spectrum. On the other hand, above the deconflnement phase transition, 
the entropy is much larger, since also quantum states in color multiplets can contribute to the 
partition function. 

In addition, flavor symmetry implies that the hadronic spectrum can be approximatively 
organized into SU(3) multiplets. Hence, since quantum states belonging to the same multiplet 
can be transformed into one another by a flavor rotation, the number of truly distinct hadronic 
wave- functions is actually much smaller than the number of hadrons recorded in the particle 
data book. Similarly, it has been experimentally observed that nearly identical native structures 
are adopted by poly-peptide chains characterized by very different amino-acid sequences. This 
is to say that the number of native structures which can be realized is much smaller than the 
number of sequences in the data bank of known proteins. 

These analogies are not accidental. They reflect the fact that the dynamics of both hadrons 
and proteins is characterized by the interplay of large correlations and quantum or thermal 
fluctuations. It is natural to address the question whether some of the powerful theoretical 
techniques developed to describe the non-perturbative internal dynamics of hadrons may be 
exported to describe the non-perturbative internal dynamics of proteins. In the next section, 
we shall see that this intuition can be made rigorous, i.e. that it is possible to formulate the 
dynamics of biomolecular systems in terms of an effective quantum many-body problem. 
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Figure 3. Comparison between the phase diagram of QCD and of proteins in solution. 



3. Path Integral Formulation of Molecular Dynamics 

The dynamics of biomolecular systems in solution is often described by means of the Langevin 
equation, which reads 

m x = -jx - WO) + f (t), (1) 

In this equation, x = (xi, X2, x/v) is a macro- vector specifying the coordinates of all the 
atomic nuclei in the molecule, U{x) is the potential energy and ^{t) is a stochastic force, which 
simulates the effects of the random Brownian collisions with the molecules in the solvent. The 
contemporary models for the potential energy U(x) — which is often improperly called the force 
field — are defined phenomenologically in order to reproduce the results of quantum chemistry 
calculations of the equilibrium configurations of small molecules [3]. For sake of simplicity and 
without loss of generality, in the following we present the formalism in the case of one particle 
in one spacial dimension. The generalization to the multi-dimensional case is straightforward. 

The acceleration term on the left hand side of Eq. introduces inertial effects which are 
damped on a time-scale t ~ m/j = rjy. For most biomolecular systems, td is of order of the 
fraction of aps, hence much smaller than the relevant microscopic time-scale associated to the 
dynamics of the torsional angles, which takes place at the ns time-scale. 

For t^> td the Langevin Eq. ([!]) reduces to 



dx 
dt 



D 



VU(x)+ri(t), 



where D 



is the so-called diffusion coefficient and 77 (t) 
that satisfy the so-called fluctuation-dissipation relationship: 



( V (t) V (t'))=2D5(t-t'). 



(2) 



is a rescaled white noise 



(3) 



It can be shown — see e.g. [4j — that the probability distribution sampled by the stochastic 
differential Eq. ^ satisfies then the Fokker-Planck (FP) equation: 



VU{x)P{x,t) +DV Z P(x,t). 



(4) 



A universal property of all the solutions of the Eq. Q is that, in the long time limit, they 
converge to the Boltzmann weight, regardless of their initial conditions. 



By performing the substitution 



P(x,t) = e ™ b t U[ - x) ip(x,t), (5) 
the FP Eq. @ can be formally re-written as a Schrodinger Equation in imaginary time: 

d 

~ *) = H eff *)> ( 6 ) 

where the effective "Quantum Hamiltonian" operator reads 

H eff = -DV 2 + V eff (x), (7) 
and V e ff(x) is an effective potential and is defined as 

V ef f(x) = [m(x)f - 2k B TV 2 U(x)] . (8) 

Hence, we have shown that the problem of studying the real time stochastic dynamics of a protein 
in solution can be mapped into the problem of determining the imaginary time evolution of an 
effective quantum many-body system. 

Based on the analogy with quantum mechanics, it is immediate to obtain a path-integral 
representation of the solution of Q, subject to the boundary conditions x(0) = Xi and x(t) = xf. 

(U(xf)-U(xi)) I I — H e fft\ 



P(x f ,t\xi) = e-^^'-^'fale-^lxf) 



= e ~^ mXf) ~ U(xi)) I "Vx(r)e Jo ^' \ 4D "—7. (9) 

JXi 

This equation expresses the fact that the Green's function of the Fokker-Planck equation is 
formally equivalent to a quantum-mechanical propagator in imaginary time. On the other hand, 
it is important to stress that the variable t which enters in these equations corresponds to the 
physical time. 



J Xi 



4. The Theory of Biological Instantons 

The path integral formulation of the stochastic dynamics of protein is particularly convenient if 
one is interested in investigating the structure of the protein folding pathways. Indeed, such a 
formalism allows to focus directly on the non-equilibrium dynamical trajectories which actually 
reach the native state Xf in a given time interval t. This way, one avoids wasting time to simulate 
trajectories which remain confined in the initial denatured state. A second important advantage 
of Eq. ^ is that it assigns a probability density to each folding pathway x(r): 

Prob.[x(r)] oc exp[—S e ff[x]]. (10) 

Clearly, the folding pathways which are most likely to be realized by a given protein are the 
solutions of the "classical" equations of motion generated by the effective action 

S ef f[x] = f Q dr ^ + Veff [x]^. (11) 

These trajectories typically overcome barriers of the effective potential V e ff(x). Hence, they 
correspond to instantons in the effective quantum theory. The so-called Dominant Reaction 



Pathways (DRP) approach developed by our group [5j EJ |9] is based on the saddle-point 
approximation of the stochastic path integral Q and focuses on such "biological instantons". 
Unlike in QCD, such an approach is not affected by infrared divergences, hence it is completely 
rigorous. The expansion parameter controlling the accuracy of this approximation is the thermal 
energy fc^T, which enters in the definition of the diffusion constant D and of the effective 
potential V e ff(x). 

In principle, the instantons can be obtained by relaxing numerically a discretized 
representation of the effective action functional S e ff[x], in complete analogy with the cooling 
procedure developed in lattice field theory. In practice, however, the simultaneous presence of 
fast and slow internal time scales (which range from fs to ns) makes such a task very challenging. 
Indeed, the discretized elementary time interval At has to be chosen much smaller than the 
smallest time scale. Hence, computing a single dominant pathway would require to find a 
minimum of a function with a gigantic number of degrees of freedom: 37V x N t , where N is the 
number of atoms in the protein — typically from several hundreds to many thousands — and Nt 
is the number of time discretization steps — typically 10 6 — . 

Fortunately, a major numerical simplification of this problem can be achieved by exploiting 



the fact that the dynamical system defined by the effective action (11) is simplectic, i.e. it 
conserves the "effective energy" Q 

E *ff = w ±2{t) - Veff[x{t)] - (12) 

Hence, rather than minimizing directly the effective action S e ff[x], it is possible to obtain the 
biological instantons using the Hamilton-Jacobi (HJ) formulation of classical mechanics. In other 
words, the trajectories obeying the classical equations of motion and subject to the boundary 
conditions x(t) = xj, and x(0) = X{ are those which minimize the effective HJ functional 

Shj[x(1)\ = -±= £ dlyjE eff (t) + V eff [x(l)}, (13) 

where dl — Vdx 2 is the measure of the distance covered by the system in configuration space, 
during the transition. Hence, the dominant reaction pathways can be viewed as the geodesic in 
a space with a curved metric. 

The major advantage of the HJ formulation is that it allows to remove the time as the 
independent variable and replacing it with the curvilinear abscissa /, which has the dimension of 
a length scale. Since there is no gap in the length scales of molecular systems, the discretization of 
the HJ is expected to converge extremely much faster than the time discretization of the effective 
action S e ff[x]. Indeed, typically 100 slices are sufficient to achieve an accurate representation 
of the path. 

The effective energy E e ff is an external parameter which determines the time at which 
each configuration of a dominant reaction pathway x(r) is visited, according to the usual HJ 
relationship 

t(x)= [ X dl 1 =. (14) 

J* J4D(E eff + V eff [x(l)]) 



Typically, one is interested in studying transitions which terminate close the local minima of 
the potential energy U(x). The residence time in such end-point configurations must be much 



longer than that in the configurations visited during the transition. From Eq. (14) it follows 



1 Note that E e ff has the dimension of a rate, therefore this quantity does not have the physical interpretation 
of a mechanical energy 



that these conditions are verified if E e fj ~ —V e ff(x ), where x Q is a configuration in the vicinity 
of a local minimum of U(x). 

In practice, computing the dominant reaction pathway connecting two given configurations 
Xi to Xf amounts to minimizing a discretized version of the effective HJ functional: 



N s -1 



71=1 



S d Hj[x(l)} = E fieff + V *ff AZ «.»+1> ( 15 ) 



where 7V S is the number of path discretization slices. Once such a path has been determined, 
one can reconstruct the time at which each of the configurations is visited during the transition, 



using Eq. (14). In particular, the time interval between the n-th and the (n + l)-th slice is 



AWl,n = , Aln+1 ' U (16) 

y/AD(E eff + V eff [x(n)]) 

where Al n +i^ n = y/(x(n + 1) - x(n)) 2 . 

4-1. Thermal Fluctuations Around the Dominant Folding Pathways 

The biological instantons encode information about the reactive dynamics, hence can be used as 
starting point to compute the folding rate, whose inverse gives the typical time a protein takes 
to reach the native state, starting from a denatured state. This is a fundamental observable in 
protein kinetics, which is accessible from fluorescence experiments. 

It turns out that, in order to obtain an equation for the rate, it is necessary to go beyond 
the lowest-order saddle-point approximation, and estimate the contribution of small thermal 
fluctuations around the dominant paths [8j. 

This corresponds to evaluating the instanton weight at one-loop level, by functionally 
expanding the effective action to quadratic order, around the instanton solution x(r): 

S eff [x] = S eff [x] + \ jf dr> jf d^ ^^^M_ {xi{T >) _ ^ (r ')) Mt) _ ^( r) ) + . . . 

" Seffm + \[ dT jl^W)-*^)) P tk'[x\ Mr)-x k (r)), (17) 
where we have introduced the fluctuation operator F[x], defined as 

S(r - r'). (18) 



1 d 2 

~2D 5 ik J^2+ 9 i d k V e ff[x(t)} 



Sxi(T')Sx k (r") 

The formal one-loop expression for the conditional probability is therefore 

-l(U(x f )-U(x t )) 

P(x f ,t\xi) ~ M = — e SefM, (19) 

ydetF[x] 

Such an expression is clearly divergent and needs to be regularized. To this end, one can 
multiply and divide by the conditional probability density of a reference system P reg .{xf,t\xi). 
For example, one can use as regulator the conditional probability associated to the diffusion in 
an external harmonic potential, for which the propagator P re g.( x o, t\%o) is analytically known. 
The result is [8]: 



PDRp{x f ,t\Xi) = Pre 9 .Mx ) e'^ U ^ U ^ e ^re 9 .]-S efM 



1 



\ det (P^l\x reg ] F\ 



x 



(20) 



where F reg is the fluctuation operator for the reference system used in the regularization. 

In practice, the determinant det (fi r eg. [x reg ] has to be evaluated numerically from 

the discretized representation of the fluctuation operators. To obtain such a representation, one 
needs to express the time derivative using discretized time intervals. It is most convenient to 



use the intervals At^+i evaluated from the dominant trajectory, according to Eq. (16). The 
fluctuation operator reads 



F\fY/ J - ~ l t D S- 



#/c,ra+l s ( 1 . 1 \ . $k,m-l 
^k,m ^T7 r ^rr + 



L -^ L m-\-l,m \ L -^ L m-\-l,m ^^m,m—l J ^^m^m—l 



d 2 V eff (x(k)) 



+ J I S k,m, (21) 

OXiOXj 

where the indexes A:, m = 1, . . . , iV s run over the path frames, while the indexes z, j = 1, . . . , d 
label the degrees of freedom of the system. 

4-2. Improved Effective Actions 

The main numerical advantage of the DRP formalism arises from the possibility of removing the 
time as a dynamical variable, and replace it with the curvilinear abscissa /. Since there is no gap 
in the characteristic length scales of molecular systems, the convergence of the discretization 
of / is usually very fast compare to that of the discretization of t. As a result, using such 
a formulation, it is possible to gain information about the reaction mechanism at a very low 
computational cost. 

On the other hand, in order to obtain information about the dynamics, one needs to compute 



the times at which each configuration is visited along the dominant path, using Eq.s (14) and 



(16). The calculation of the time intervals At^+i is also needed in the discretized representation 



of the fluctuation operator (21), which enters in the calculation of the order fc^T corrections 
arising from non-equilibrium stochastic fluctuations around the dominant path. 

Clearly, the trick to remove time by switching to the HJ formulation does not help in the 
evaluation of explicitly time-dependent observables. As a result, in order to achieve an accurate 
description of the dynamics, or in order properly take into account of the effects of fluctuations, 
one would need to use a large number of path frames, with a consequent significant increase of 
the computational cost of a DRP simulation. 

Fortunately, the convergence of the calculation of the time intervals in the DRP approach can 
be greatly improved by adopting the effective stochastic theory (EST) developed in [10J. The 
main idea of the EST is to exploit the gap in the internal time scales in order to analytically 
perform the integral over the fast Fourier components of the paths x(r) which contribute to the 
path integral ([9]), using Wilson's approach. Through such a procedure, the effects of the fast 
dynamics is rigorously and systematically integrated out by renormalizing the effective potential: 

V eff (x) -+ V e f f T (x) = V eff (x) + V e I } f (x). (22) 

where, to leading order in the effective theory one has 

y3M = DAt ;X~ b) v2v ^ x)+ --- {23) 

In such an Eq., At c is a cut-off time scale which must be chosen much smaller than the fastest 
internal dynamical time scale and b is a parameter which defines the interval of Fourier modes 
which are being analytically integrated out. Typically, for molecular systems At c ~ 10 _3 ps and 



b rsj 10 [TT]. The dots in Eq. (23) denote higher order correction in an expansion in the ratio 



of slow and fast time scales (so-called slow-mode perturbation theory) and can be found in the 
original publication [TO] . 

The EST generates by construction the same long-time dynamics of the original "bare" 
theory, but has a lower time resolution. In the context of MD simulations, this implies that the 
EST can be integrated using much larger discretization time steps [IT] . In the context of the 
DRP simulations, the utility of the EST resides in the fact that much fewer path discretization 
time steps are required in order to achieve a convergent calculation of the time interval from the 
dominant path, through Eq. (14). 



4-3. Generalization to Ab-initio Quantum Chemistry Calculations 

So far we have discussed the DRP formalism in the specific context of protein folding. On the 
other hand, this theory applies in general to any physical system described by the over-damped 
Langevin Eq. ([2]). In particular, it can be applied to investigate rare chemical reactions taking 
place in solution [12J. Typically, such reactions involve large modifications of the electronic 
structure, e.g. associated to the breaking or forming of covalent bonds. Hence, in investigating 
such processes, a classical description based on some phenomenological expression for the 
molecular potential energy U{x) entering in the Langevin equation is no longer appropriate, 
and one must adopt a sub-atomic approach, in which the quantum electronic degrees of freedom 
are explicitly taken into account. 

In the Born-Oppenheimer approximation [3], this is done in two steps. First, the Schrdinger 
equation for the electronic degrees of freedom is (approximatively) solved holding fixed all the 
nuclear coordinates. The corresponding ground-state energy is then interpreted as the potential 
molecular energy U(x) and is used to derive the effective potential V e ff(x), which enters the 
DRP formulas. This approach allows to rigorously take into account of the quantum effects on 
the electron dynamics, to all order in % and to leading order in the ratio m/M between the mass 
of the electron m and the typical mass of the atomic nuclei M. 



4-4- Quantum Corrections to the Diffusive Dynamics of the Atomic Nuclei 
In the discussion performed so far, the dynamics of the atomic nuclei has been assumed to be 
entirely classical, and to evolve according to the Langevin Eq. Q. Such an approximation 
is certainly reliable for most atomic species which are found in biomolecular systems, such as 
carbon or oxygen. On the other hand, it becomes questionable for the lightest species, notably 
hydrogen. Hence, in the study of reactions involving a large number of such atoms, the quantum 
corrections to the stochastic dynamics of nuclei should be consistently taken into account. 

To tackle this problem, we have derived the multidimensional generalization [15] of the 
semiclassical extension of the Langevin Eq. ^ and of the corresponding FP equation Q, in 
which quantum corrections are systematically included to order f? [131 (HI- From these equations 
it is immediate to compute the quantum corrections to the DRP equations [15] . We found that 
the effective action receives contribution from an additional term in the form (adopting again 
for simplicity a one-dimensional notation): 

S? ff [x] = J* dr V e Q ff [x(r)}, (24) 

where 

V?„(x) ~ ^J^\VU(x)\ 2 V 2 U(x) (25) 

is the correction to the effective potential and A = 12 k B Tm- ls a cnarac teristic parameter 
which determines the size of quantum fluctuations. In conformational transitions driven by 
the formation of hydrogen bonds, the quantum corrections to the dominant reaction pathways 
have been found to give raise to sizable effects, see Ref. [15] • 
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Figure 4. Evolution of the end-point distance of the protein during a folding reaction, averaged 
over folding trajectories obtained in DRP and MD simulations, respectively. The dashed line 
represents a typical trajectory. 

5. Some Applications 

The DRP approach has been extensively tested and validated on toy systems [8j and 
simplified protein models [HUE], and then applied to investigate realistic transitions, including 
conformational transitions of peptides [HI [18] , chemical reactions [12] and realistic protein folding 
reactions [2T] . In this section, we briefly review some of these studies. 

5.1. Validation of the DRP Theory on a Coarse- Grained Model 

The most straightforward way to assess the accuracy of the DRP approach consists in comparing 
its predictions with those obtained directly from MD simulations. Unfortunately, for realistic 
models MD simulations of protein folding reactions are not presently feasible. 

Hence, we have validated the DRP approach using much simpler coarse-grained models, with 
protein-like properties. For example, the model used in Ref. [17J displays a folding transition 
with a well-defined and unique native state. The predictions for the time evolution of the 
distance between the end-points of the chain (a typical observable in protein folding) obtained 
in the DRP and MD approaches is compared in Fig. [4j We see that the biological instanton 
standpoint provides a realistic description of the reactive dynamics. 

5.2. Atomistically Detailed Protein Folding Simulations 

The computational efficiency of the DRP approach allows to investigate a large class of reactions 
which cannot be studied on the existing supercomputer using standard MD simulations. Only 
very few examples of successful MD simulations of protein folding reactions have been reported 
to date. These studies were either performed on a special purpose supercomputing machine in 
which all the key MD algorithms are implemented at the hardware level [19] or through very long 
large-scale simulations performed on several thousands of CPU's, distributed world wide [20J. 

On the other hand, using the DRP method [21] it was possible to simulate the same reaction 
at a comparable level of statistics in less than a single day on only 32 CPU's. The DRP results 
were found to be consistent with those reported in [15] . and provided an explanation of the 
discrepancy between the results of Ref. [19J and Ref. [20]. An illustration of one of the folding 
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Figure 5. Folding of a peptide obtained from quantum DRP calculations. The shaded area 
represents the distribution of electric charge. 

trajectories obtained in the DRP approach is displayed in Fig. [I] Simulations on much larger 
systems (unaccessible to MD simulations) are now being performed. 

5.3. Folding of a Small Peptide from Ab-initio Quantum Mechanical Simulations 
Using the DRP approach we have performed the first study of the folding of a peptide 
based on ab-initio Quantum Mechanical calculations, using the Born-Oppenheimer approach 
described in section [43} This type of analysis could never performed using MD or Car-Parinello 
approaches [3j. The results of our simulations are schematically represented in Fig. [5j This 
study allowed for the first time to access the reliability of the existing classical force fields in 
predicting non-equilibrium reactive trajectories. This is an important and non-trivial issue, 
since the phenomenological force fields are fitted in order to reproduce quantum mechanical 
calculations on static configurations, or in equilibrium condition. Hence, there is in principle no 
guarantee that they should be accurate in also in the regions of phase-space which are mostly 
visited by the non-equilibrium reactive trajectories. On the other hand, our studies have revealed 
that classical calculations give results that are indeed in good agreement with those of quantum 
ab-initio simulations. This study also allowed to quantify to what extent the electronic structure 
of the chain is modified during the reaction. 

6. Conclusions 

Nuclear physics provides an excellent training ground where developing advanced theoretical 
approaches to interacting many-particle systems. Indeed, a number of schemes and 
approximations originally proposed in this context [22] have later become standard tools in 
the field of strongly correlated electrons, quantum liquids, cold atoms, .... 

In particular, in this talk we have reviewed our recent attempt to export functional integral 
techniques originally derived in QCD to investigate the non-equilibrium diffusive dynamics of 
biomolecules such as proteins. We have argued that such an approach opens the door to the 
study of a large number of processes which simply could not be investigated using standard MD 
simulations. 

Interdisciplinary studies of this type may also emphasize the importance of pushing forward 
the fundamental research in nuclear theory. Indeed, they show that this activity is valuable also 
in a perspective which goes beyond the traditional perimeter of nuclear physics. 
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